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The lateral distribution function of high energy muons in EAS around the knee (5.9 < lgiV e < 7.1) has been 
measured for near vertical showers (9 < 20° , effective muon threshold energy is 230 GeV) . The measurements have 
been performed at the Baksan Underground Scintillation Telescope (BUST). The electromagnetic component is 
measured by the "Andyrchy" EAS array, located above the BUST. The knee in EAS size spectrum is found to 
be at lg N e ~ 6.3. The experimental results are compared with Monte Carlo simulations. 



1. Introduction 

Measuring the lateral distribution function of 
high energy muons and its dependence on the 
shower size provides additional information for 
studying the mass composition of the primary 
cosmic radiation as well as for choosing a model of 
high energy hadron interactions. In this work we 
present results of measurements of the muon lat- 
eral distribution functions in EAS for the follow- 
ing three regions around the knee: 5.9 < lg N e < 
6.3, 6.3 <lgN e < 6.5 and 6.5 < lgN e < 7.1 (the 
muon energy threshold is 230 GeV) . The measure- 
ments have been performed using the facilities of 
the Baksan Neutrino Observatory: the BUST and 
the "Andyrchy" EAS array Q]2]. The live data 
taking time is 24460.8 hours. The knee position 
in the EAS spectrum corresponding to the con- 
sidered zenith angles 8 < f0° is found to be at 
lgiV e 6.3. We have also simulated the develop- 
ment of EAS in the Earth atmosphere by means 
of the CORSfKA package (version 6720). As the 
high and low energy hadronic interaction mod- 
els we used QGSJET 01C and GHEfSHA 2002d, 
respectively [3]. The simulations have been car- 
ried out for primary protons and iron nuclei and 
the obtained muon lateral distribution function is 
compared with the experimental one. 
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2. Experiment 

The "Andyrchy" EAS array is located on the 
slope of the mountain Andyrchy above the BUST 
and consists of 37 plastic scintillation detectors of 
1 m 2 area each. The distance between the detec- 
tors is about 40 m in projection to the horizon- 
tal plane and the overall area of the installation 
is 4.5 • 10 4 m 2 . The detectors are arranged in 
such a way that the central one is just above the 
BUST at a vertical distance of about 360 m which 
corresponds to 2060 m above the sea level. The 
shower trigger condition requires four detectors 
to fire within 3 microseconds. The trigger's rate 
is about 9 s _1 . The array and its characteristics 
are described in more details in g]. 

The BUST is a four-floor building with 
16.1 x 16.1 x 11.2 dimensions located in a mine 
at the effective depth 850 hg/cm 2 . The floors as 
well as the four vertical sides of the building are 
fully covered with 3f 50 liquid scintillation detec- 
tors. Each of the detectors has the dimensions 
(0.7 x 0.7 x 0.3) m 3 . The telescope allows to 
determine the number of the passing muons (1 
- 200), their coordinates (with 0.7 m accuracy) 
and the arrival direction (with 1.5 degree accu- 
racy). The coincidence trigger rate of the BUST 
and "Andyrchy" is about 0.1 s" 1 . 

The size, axis position and the EAS arrival di- 



1 



2 



rection are determined using the "Andyrchy" ar- 
ray data; the BUST data are used to determine 
the number of muons crossing it. In the present 
analysis, only the showers with axes in central 
part of the installation (the distance from the 
center is not larger than 50 m) and with 5.9 < 
lgiV e < 7.1 were taken into account. The accu- 
racy of the determination of N e for such showers 
is not worse than 12%. For more details on iden- 
tification of EAS parameters at the "Andyrchy" 
detector array see, for example, [51617] . 

3. The fraction of muons in the telescope 

The underground telescope measures only a 
part of the total number of muons in EAS and 
the uncertainty in the determination of the EAS 
axis position at the observation level is compara- 
ble with the size of the BUST. In our experiment 
the fraction of muons in the telescope A(R) was 
measured as a function of the distance R between 
the center of the telescope and the EAS axis for a 
set of showers with a given value of N e [819] . This 
fraction for a given range of N e is defined by 
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Figure 1. The mean number of muons in the tele- 
scope as a function of the distance R between 
the center of the telescope and the shower axis 
for showers with 5.9 < lgA^ < 6.3. The exper- 
imental measurements are depicted by the dots. 
The curves correspond to the calculations with 
r = 4.6 m using equation [5] and for (Tr — m 
(curve 1) and <tr = 12 m (curve 2). 



where n(R) is the mean number of muons in the 
telescope at the distance R, N fi (N e ) is the mean 
number of muons in EAS. The mean number of 
muons in the BUST is determined as follows. 
Events within a given range of N e are grouped 
according to the distance to BUST center with 
step SR = 10 m. For each group, the number of 
muons is 
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where Ki — K{Ri) is the number of EAS in the 
i-th group, and my is the total number of muons 
in the BUST in the i-th group for the j-th EAS. 

Thus, the mean number of muons for a given 
range of N e can be written as 
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The mean number of muons in a shower is then 
calculated as 



N, 



n{Ri) -Sr(Ri), 
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where St — 200 m 2 is the effective area of the 
telescope, S r (Ri) is the area of the ring of radius 
Ri with SR = 10 m. 

We have mentioned above that the fraction of 
muons with > 230 GeV in the telescope as 
a function of R was obtained for the following 
three energy ranges about the knee in the EAS 
spectrum. For the first one, before the knee 
5.9 <\gN e < 6.3 (\gN e = 6.1), the mean num- 
ber of muons in a shower is N ^ = 100 ± 2; for 
6.3 < lgA^ e < 6.5 (lgiV e = 6.4), thisjs JV P = 
184 ± 9; and for 6.5 < lgAT e < 7.1 (\ g N e = 6.7) 
= 334 ± 19. The experimentally measured 
dependence A(R) for 5.9 < lgAT e < 6.3 is shown 
in Fig. [Has black dots. 
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4. The muon lateral distribution function 



The fraction of muons in the telescope A(R) 
depends only on the muon lateral distribution 
function f(r) and the geometry of the installa- 
tion. For distances much larger than the typical 
size of the installation (R 3> \/St), the fraction 
of muons is A(i?) « St ■ f{R)- In the analysis of 
the experimental data, the following form of the 
function was used 



a + b- 



(5) 



which, according to the BUST data on the muon 
groups, sufficiently reproduces the lateral distri- 
bution of high energy muons with the threshold 
energy-dependent parameter r , when a = 0.58, 
b = 0.64, c = 0.47 [lOlllj . The same function fits 
the lateral distribution of muons with the thresh- 
old energy 230 GeV in the simulated showers as 
well. However, in order to keep the good agree- 
ment at the primary energies above 10 5 GeV per 
nucleon, one has to change the values of the pa- 
rameters to: a — 0.65, b — 0.72, c = 0.4, while ro 
depends only on the primary energy per nucleon. 

The calculation of A(i?) has been performed 
by the Monte Carlo method for a few values of 
r with taking into account the uncertainty in 
determination of the shower axis which, at the 
BUST level, is <tr ~ 12 m (the main contribu- 
tion is from the determination of the EAS arrival 
direction due to the large distance between the 
installations). 

In Fig. [T]we show the experimentally measured 
A(i?) for 5.9 < lgN e < 6.3 in comparison with 
the calculations at ro = 4.6 m for two values of 
<jr: m (no uncertainty in the determination of 
the arrival direction) and 12 m. This comparison 
gives the value of ro for each of the ranges of N e 
(Fig. 02): for lgJVe = 6.1 r = 4.6 ± 0J5 m, for 
\gN e = 6.4 r = 4.0 ± 0.8 m and for lgN e = 6.7 
r = 3.9 ±0.5 m. 

The dependence of the experimental values of 
ro on lgiV e in comparison with the correspond- 
ing calculations for protons and iron nuclei is dis- 
played in Fig. O In this case the mean values of 
N e are taken as the means for every fixed primary 
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Figure 2. Dependence of the parameter ro in the 
muon lateral distribution function on the EAS 
size. Points - experiment, lines - simulations for: 
1 - primary protons, 2 - primary iron nuclei. 



energy, without taking fluctuations and the form 
of the primary spectrum into account. 

This work was supported by the "Neutrino 
Physics" Program for Basic Research of the Pre- 
sidium of the Russian Academy of Sciences and 
by "State Program for Support of Leading Sci- 
entific Schools" (project no. NSh-321.2008.2). 
This work was also supported in part by the Rus- 
sian Foundation for Basic Research (grants 06-02- 
16135 and 08-07-90400). 

REFERENCES 

1. E.N. Alexeyev et al., Izv. RAS, Ser. Phys. 57 
(4) (1993) 99. 

2. E.N. Alexeyev et al., Proc. 16th ICRC, Kyoto, 
1979, v.10, p.276. 

3. D. Heck et al., Report FZKA 6019 (1998), 
Forschungszentrum, Karlsruhe. 

4. V.B. Petkov et al., Instruments and Experi- 
mental Techniques, 2006, v. 49, No. 6, p. 785. 

5. V.B. Petkov et al., Izv. RAS, Ser. Phys. 66 
(11) (2002) 1560. 



4 



6. V.B. Petkov et al., Proc. 28th ICRC, 
Tsukuba, 2003, v.l, p.65. 

7. A.V. Voevodsky et al., Izv. RAS, Ser. Phys. 
61 (3) (1997) 496. 

8. A.E. Chudakov et al., Proc.25 ICRC, Durban, 
1997, v.6, p.173. 

9. A.E. Chudakov et al., Proc. Xth International 
School "Particles and Cosmology", Moscow, 
2000, p.303. 

10. V.N. Bakatanov et al., Proc. 18 ICRC, Ban- 
galore, 1983, v.ll, p.453. 

11. A.V. Voevodsky et al., Yad. Fiz. 56 (1993) 
143. 



